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Abstract 

This book chapter introduces to the concept of weak chaos, aspects of its ergodic theory 
description, and properties of the anomalous dynamics associated with it. In the first half 
of the chapter we study simple one-dimensional deterministic maps, in the second half basic 
stochastic models and eventually an experiment. We start by reminding the reader of fun¬ 
damental chaos quantities and their relation to each other, exemplified by the paradigmatic 
Bernoulli shift. Using the intermittent Pomeau-Manneville map the problem of weak chaos 
and infinite ergodic theory is outlined, defining a very recent mathematical field of research. 
Considering a spatially extended version of the Pomeau-Manneville map leads us to the 
phenomenon of anomalous diffusion. This problem will be discussed by applying stochas¬ 
tic continuous time random walk theory and by deriving a fractional diffusion equation. 
Another important topic within modern nonequilibrium statistical physics are fluctuation 
relations, which we investigate for anomalous dynamics. The chapter concludes by showing 
the importance of anomalous dynamics for understanding experimental results on biological 
cell migration. 
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1 Introduction 


Deterministic dynamical systems involving only a few variables can exhibit complexity remi¬ 
niscent of many-particle systems if the dynamics is chaotic , as is quantified by the existence 
of a positive Lyapunov exponent [Sch89b, Ott93, Bec93, A1197]. Such systems, which may 
be called small because of their small number of degrees of freedom [Ivlal2b], can display 
an intricate interplay between nonlinear microscopic dynamical properties and macroscopic 
statistical behavior leading to highly non-trivial fluctuations of physical observables. This 
becomes particularly interesting in nonequilibrium situations when these systems are ex¬ 
posed to external gradients or fields. Despite their complexity, examples of these systems 
are still amenable to detailed analysis by means of dynamical systems theory in combination 
with stochastic theory. Hence, they provide important paradigms to construct a theory of 
nonequilibrium statistical physics from first principles: Based on the chaotic hypothesis, 
which generalizes Boltzmann’s ergodic hypothesis, SRB measures were studied as nonequi¬ 
librium equivalents of the Gibbs ensembles of equilibrium statistical mechanics. This novel 
approach led to the discovery of fundamental relations characterizing nonequilibrium trans¬ 
port in terms of microscopic chaos [Dor99, Gas98, Kla07, Cas08], such as formulas expressing 
transport coefficients in terms of Lyapunov exponents and dynamical entropies, equations 
relating the nonequilibrium entropy production to the fractality of SRB measures, and fluc¬ 
tuation relations, which are now widely studied as a fundamental property of nonequilibrium 
processes [Eva02, Bus05, Kla07, Klal2b]. 

The interplay between these different levels of description in modern nonequilibrium sta¬ 
tistical mechanics is illustrated by the second column in Fig. 1.1, in analogy to the theory 
of equilibrium statistical mechanics sketched in the first column. As is represented by the 
third column, however, more recently scientists learned that random-looking evolution in 
time and space also occurs under conditions that are weaker than requiring a positive Lya¬ 
punov exponent [Aar97, ZasOl]. ft is now known that there is a wealth of systems exhibiting 
zero Lyapunov exponents, meaning that the separation of nearby trajectories is weaker than 
exponential. This class of dynamical systems is called weakly chaotic. Examples include 
maps with indifferent fixed points, polygonal particle billiards, and Hamiltonian systems 
with sticky islands in phase space [Aar97, ZasOl, Kla07, Kla08]. 

Weakly chaotic systems exhibit anomalous dynamics characterized by novel properties such 
as ageing, which reflects an extremely weak relaxation towards equilibrium involving more 
than one time scale in the decay of correlations. Other surprising properties are the existence 
of Levy-type probability distributions obeying generalized central limit theorems [Shl93, 
Kla96a] and the non-equivalence of time and ensemble averages, called weak ergodicity 
breaking [Ste09] . These physical phenomena were observed experimentally in a wide variety 
of systems, such as in the anomalous statistics of blinking quantum dots, in the anomalous 
diffusion of atoms in optical lattices, in plasma physics and even in cell and animal migration 
[MetOO, Met04, Kla08, Ste09]. 
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Figure 1.1: Conceptual foundations of a theory of nonequilibrium statistical physics based on 
dynamical systems theory by motivating the topic of this book chapter, which is represented 
by the third column. 


Recent work in ergodic theory, on the other hand, has led to mathematically rigorous results 
about some of the physically relevant phenomena mentioned above. It turned out that 
there is an intimate connection between the mechanism generating weakly chaotic dynamics 
and the existence of non-normalizable, so-called infinite invariant measures [Aar97, AarOO, 
Tha06]. The ergodic theory of generalized random walks driven by weak chaos and of 
other systems exhibiting infinite invariant measures, which is called infinite ergodic theory , 
has thus the potential of providing a sound mathematical basis for some of the physical 
phenomena displayed by anomalous dynamics. 

This book chapter gives a brief introduction to important aspects of the above topics in four 
sections: As a warm-up, the beginning of Section 2 briefly reminds us of the concept of deter¬ 
ministic chaos in simple dynamical systems as quantified by a positive Lyapunov exponent. 
On this basis, we will introduce to the phenomenon of weak chaos, and the idea of infinite 
ergodic theory will be outlined. The chapter concludes by putting different forms of deter¬ 
ministic chaos into perspective. Section 3 relates these concepts and ideas to the problem of 
anomalous diffusion in deterministic systems. Here we make a transition to stochastic theory 
by studying these systems also from a stochastic point of view. For this purpose we use a 
generalization of ordinary random walk theory, called continuous time random walk theory. 
In a scaling limit, this theory leads to generalized diffusion equations involving fractional 
derivatives. Section 4 introduces to the topic of fluctuation relations, which generalize the 
Second Law of Thermodynamics and other fundamental thermodynamic relations to small 
systems far away from equilibrium. After discussing transient fluctuation relations for a very 
basic type of stochastic dynamics as an example, we explore the validity of such relations for 
generalizations of this dynamics yielding anomalous diffusion. In section 5 we relate this line 
of theoretical reasoning about anomalous dynamics to biophysical reality by studying the 
case of biological cell migration. After briefly introducing to the problem of cell migration, 
we report experimental results on fundamental statistical physical properties of migrating 
cells, extracted from statistical data analysis. We conclude this section with a stochastic 
modeling of these experimental results by using a generalized, fractional Fokker-Planck type 
equation. We summarize our discussion of this book chapter in the final section 6. 




The title of this review is inspired by a conference that the author had the pleasure to 
organize together with R.Zweimuller, E.Barkai, and H.Kantz at the Max Planck Institute 
for the Physics of Complex Systems, Dresden, in Summer 2011, which bears exactly the 
same title [Klall]. However, naturally this chapter reflects the author’s very personal take 
on this topic, and his own research. The subsequent second section is to some extent based 
on the review Ref. [KlalO] by combining it with ideas from Refs. [I\la07, How09]. The third 
section bnilds on Refs. [Kor05, Kor07]. The fourth section incorporates material from the 
review Ref. [Ivlal2a] and from Ref. [Che09], the fifth one draws on Ref. [Die08]. 


2 Chaos and anomalous dynamics 


In this section we focus on purely deterministic dynamics modeled by two simple but 
paradigmatic one-dimensional maps: the famous Bernoulli shift, as a model for strong chaos 
characterized by a positive Lyapunov exponent, and the Pomeau-Manneville map, as an ex¬ 
ample exhibiting weak chaos with zero Lyapunov exponent. We start by briefly reminding 
of basic concepts of dynamical systems theory and ergodic theory such as Lyapunov expo¬ 
nents, ergodicity, SRB measures, and Pesin’s theorem, illustrated for the Bernoulli shift. 
Ref. [KlalO] provides a more tutorial exposition of most of these ideas. By switching to the 
Pomeau-Manneville map we find that generalizations of these concepts are needed in order 
to describe the model’s weakly chaotic dynamics. This motivates the mathematical problem 
of infinite ergodic theory, which is intimately related to defining suitably generalized chaos 
quantities assessing weak chaos, and a generalization of Pesin’s theorem. In the final part of 
this chapter we propose a generalized hierarchy of chaos, based on the existence of different 
types of stretching between two nearby trajectories, which we use to characterize chaotic 
dynamics. 


2.1 Deterministic chaos in a simple map 

The main vehicle of our approach in this and the next section are one-dimensional time- 
discrete maps F : J —>■ J , JCi obeying 

x n+1 = F(x n ) , n G N 0 , (2.1) 

which defines the equations of motion of our deterministic dynamical systems. For a given 
initial condition xo we have x n = F n (x o). A particularly simple example of F are piecewise 
linear maps, such as the paradigmatic Bernoulli shift [Sch89b, Ott93, A1197, Dor99] 

B : [0,1) [0,1) , B(x) := 2x mod 1 = { 3 * <1 < j (2.2) 

depicted in Fig. 2.1. This simple system exhibits a very complicated dynamics governed by 
sensitivity to initial conditions, as can be quantified by calculating its Lyapunov exponent 
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Figure 2.1: The Bernoulli shift Eq. (2.2) and two trajectories starting from two nearby 
initial conditions xq and x' 0 = xq + Axo displaced by Axq -C 1. 

[Ott93, Rob95]: Consider two points that are initially displaced from each other by Ax 0 : = 
\x' Q — Xq\ with Axq “infinitesimally small” such that Xq,x' q do not hit different branches of 
the Bernoulli shift B{x) around x = 1/2 . 1 We then have 

Ax n := \x' n — x n \ = 2Ax n -\ = 2 2 Ax n _2 = ... = 2 n Axo = e nlu2 Axq . (2.3) 

We thus see that there is an exponential separation between two nearby points as we follow 
their trajectories, where the rate of separation A(xo) := In 2 is the (local) Lyapunov exponent 
of B{x). Since A(xo) > 0 , the system displays an exponential dynamical instability and is 
hence called chaotic (in the sense of Lyapunov) [Rob95, Ott93, A1197, Bec93]. 

Writing down the analogue of Eq. (2.3) for a given differentiable map F, we get 

Ax n = K-x n \ = | F n (x' 0 ) - F n {x 0 )| =: e^Axo (Ax 0 A 0) , (2.4) 

which we can take as the definition of the Lyapunov exponent A(xo) that comes in as the 
exponential stretching rate on the right hand side. Solving this equation for A(xo) by using 
the chain rule, it is not too hard to see [A1197] that this simple procedure of calculating A 
can be generalized in terms of the time (or Birkhoff) average 

1 n— 1 

A(x) = lim — ln|F , (x;)| (2.5) 

n—>oo Ti * J 
»=0 

with x = Xq. If the dynamical system defined by the map F is ergodic, the time average 
does not depend on the initial condition for a typical x, X = X(x) = const. It can be 
shown that the Bernoulli shift is ergodic [Dor99], and indeed, following Eq. (2.5), for B 
we trivially find that A = In 2 for all x. In particular, according to Birkhoff’s theorem 
[Dor99, Kat95, Arn 68 , Tod92], for ergodic systems the time average is equal to the ensemble 
average, which for the Lyapunov exponent of one-dimensional maps reads 

A = (In \F'(x)\)fj,* := j d^*hi\F\x)\ . (2.6) 


1 This condition could be eliminated by defining a metric on a circle [A1197] . 
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Here /i* is the invariant measure of the map. If the map exhibits an SRB measure [You02, 
Eck85, BalOO], we have 

dn* — p*(x) dx , (2.7) 

where p*(x) holds for the invariant density of the map. That is, the measure /T has the nice 
property that it can be obtained by integrating a density, 


p*(A) = dx p*(x) , A C J , 

J A 


( 2 . 8 ) 


which simplifies the calculation of the ensemble average Eq. (2.6). For the Bernoulli shift it 
is not too difficult to see [Bec93] that, for typical initial conditions, the invariant density is 
p*(x) = 1. By combining Eqs. (2.6) and (2.7), we get 


A 


dxp* ( x ) In 2 = In 2 


(2.9) 


This result is equal to the time average calculated above and confirms the result obtained 
from our handwaving argument Eq. (2.3). 

Lyapunov exponents are not the only quantities assessing the chaotic character of a dy¬ 
namical system. Pesin’s Theorem [Dor99, Eck85, You02] states that for closed C 2 Anosov 
[Eck85, Dor99] systems the Kolmogorov-Sinai (or metric) entropy h^s is equal to the sum of 
positive Lyapunov exponents. For one-dimensional maps that are expanding [A1197, Bec93], 


VxeJ |F'(x)l>l, 


( 2 . 10 ) 


this theorem boils down to 
where [Ott93, Dor99] 

hxs 


A — hxs , 


lim- 

n—>oo Tl 


Y p*(w)hip*(w). 
we{wp} 


( 2 . 11 ) 

( 2 . 12 ) 


Here p*(w ) is the SRB measure of an element w of the partition {H 7 / 1 }, and n defines 
the level of refinement of the partition. Note that in Eq. (2.12) we have assumed that the 
partition is generating [Eck85, Kat95, Bad97]. If hxs > 0 one sometimes speaks of measure- 
theoretic chaos [Bec93]. For the Bernoulli shift it is not too hard to calculate h^s from first 
principles leading to h KS = In 2 [Ott93, Bec93], which combined with our previous result 
for the Lyapunov exponent is in line with Pesin’s theorem. This theorem can be formulated 
under weaker assumptions, and it is believed to hold for a wider class of dynamical systems 
than stated above. We remark that typically the KS-entropy is much harder to calculate 
for a given dynamical system than Lyapunov exponents. Hence, Pesin’s theorem is often 
employed in the literature for indirectly calculating the KS-entropy. 


2.2 Weak chaos and infinite ergodic theory 

Let us now consider a nonlinear generalization of our previous piecewise linear model, which 
is known as the Pomeau-Manneville map [Pom80], 


P a ,z{x) = x + ax z mod 1 


(2.13) 
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P ( 



X 


Figure 2.2: The Pomeau-Manneville map Eq. (2.13) for a = 1 and z = 3. Note that there 
is a marginal fixed point at x = 0 leading to the intermittent behavior depicted in Fig. 2.3. 

see Fig. 2.2, where following Eq. (2.1) the dynamics is defined by x n+ \ = P a ,z{x n ). This 
map has the two control parameters a > 1 and the exponent of nonlinearity z > 1. For 


a = 1 and z — 1 the map reduces to the Bernoulli shift Eq. (2.2), for z > 1 it provides a 
nontrivial nonlinear generalization of it. The nontriviality is due to the fact that in this case 
the stability of the fixed point at x = 0 becomes marginal (sometimes also called indifferent, 


or neutral), P' az fiS) = 1. This implies that the map is non-hyperbolic, because [Dev89], 


flN > 0 such that VxVn > N \{Pa, z )'i x )\ ^ 1 > 


(2.14) 


which is related to the fact that the map is not expanding anymore according to Eq. (2.10). 
Since the map is smooth around x — 0, the dynamics resulting from the left branch of the 
map is determined by the stability of this fixed point, whereas the right branch is still of 
‘Bernoulli shift-type’ generating ordinary chaotic dynamics. There is thus a competition in 
the dynamics between these two different branches as illustrated by the time series displayed 
in Fig. 2.3: One can observe that long periodic ‘laminar phases’ determined by the marginal 
fixed point around x = 0 are interrupted by ‘chaotic bursts’ reflecting the Bernoulli shift- 
type part of the map with slope a > 1 around x — 1. This phenomenology is the hallmark 
of what is called intermittency [Sch89b, Ott93] . 

This seemingly small nonlinear modification of the Bernoulli shift has dramatic consequences 
for the whole dynamics of the new system. We discuss them step by step following our 
exposition of the Bernoulli shift dynamics in Section 2.1: The invariant density of the 
Pomeau-Manneville map can be calculated to [Tha83, Zwe98, Kor09, KorlO, How09] 


p*(x ) ~ x 1 z (x —)■ 0) . 


(2.15) 


This singularity reflects the stickiness of trajectories to the marginally unstable fixed point 
at x = 0. Correspondingly, the measure obtained by integrating this density is non- 
normalizable for z > 2 yielding the infinite invariant measure 



dy p*(y) —> oo (x —> 0) . 


i 


(2.16) 
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Figure 2.3: Phenomenology of intermittency in the Pomeau-Mannevillc map Fig. 2.2: The 
plot shows the time series of position x n versus discrete time step n for an orbit generated 
by the map Eq. (2.13), which starts at a typical initial condition xq. 

The branch of erogdic theory exploring the ergodic properties of infinite measure preserving 
dynamical systems is thus called infinite ergodic theory ; see Refs. [ThaOl, AarOO, Zwe09] 
for introductions to this topic and [Aar97] for an in-depth mathematical treatment. The 
marginal fixed point has also an impact on the dispersion of nearby trajectories, which can 
be calculated to [Gas88, Kor09, KorlO, How09] 

Ax n ~ exp Ax 0 (z > 2 ) . (2-17) 

In contrast to the Bernoulli shift, which according to Eq. (2.3) exhibits exponential sensi¬ 
tivity to intial conditions, here we thus have a weaker stretched exponential sensitivity. By 
repeating the calculation leading to Eq. (2.5), it is not hard to see that Eq. (2.17) yields a 
zero Lyapunov exponent, 

A = 0, (2.18) 

despite the fact that Fig. 2.3 displays irregular dynamics. Dynamical systems where the 
separation of nearby trajectories grows weaker than exponential, which implies that the cor¬ 
responding Lyapunov exponents are zero, have been coined weakly chaotic [ZasOl, Gal03, 
vB04, Art05]. We remark, however, that this denotation is not used unambiguously in the 
literature. Most importantly, the standard definitions of Lyapunov exponents for expand¬ 
ing and hyperbolic systems Eqs. (2.5),(2. 6 ) yield no good indicators of irregular dynamics 
anymore, because they do not capture the sub-exponential dispersion of trajectories. It is 
thus desirable to come up with generalized definitions of ordinary chaos quantities, which 
enable us to still assess this different type of chaotic behavior by calculating quantities that 
yield non-zero values. 

The way to achieve this goal is shown by advanced concepts of infinite ergodic theory and 
corresponding generalized ergodic theorems. Recall that Birkhoff’s theorem implies that for 
observables which are Lebesgue integrable, / e L 1 , we have [Dor99, Kat95, Arn 68 , Tod92] 

-j 71—1 

-5^/0d) = </v • 

i=0 


(2.19) 
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However, it turns out that for z > 2 the Birkhoff sum on the left hand side does not 
converge anymore. Surprisingly, it becomes a random variable that depends on initial 
conditions, and the equation breaks down. This non-equivalence between time and ensem¬ 
ble averages became known as weak ergodicity breaking in the physics literature, see, e.g., 
Refs. [Reb08, Kla08, Ste09] and further references therein. It was observed experimentally 
in the anomalous statistics of blinking quantum dots and plays also a crucial role for the 
anomalous diffusion of atoms in optical lattices [Kla08, Reb08, Ste09] . Note that physicists 
typically refer to ergodicity as the equality between time and ensemble average, whereas 
mathematicians usually define ergodicity via indecomposability [Arn 68 , Kat95]. Eq. (2.19) 
then follows from this definition by using Birkhoff’s theorem. This should be kept in mind 
when referring to a weak ergodicity breaking. 

In case of z > 2 and / e L 1 for our map, the nature of the breakdown of Eq. (2.19) is 
elucidated by the Aarons on-Darling-Kac theorem [ZweOO, Tha06, Zwe09] 


1 

d n 


n— 1 

E f(Xi) Ma(f)fj* (n -t oo) , 

i =0 


( 2 . 20 ) 


where the arrow holds for convergence in distribution. Here A4 a , a G [0,1], denotes a 
non-negative real random variable distributed according to the normalized Mittag-Leffler 
distribution of order a, which is characterized by its moments 


(Ml) 


, (r(i + q)) r 

T(1 + rot) 


r > 0 . 


( 2 . 21 ) 


For the Pomeau-Manneville map P a , z one can prove [ZweOO] that a n ~ n a with a := l/(z— 1 ). 
Integrating Eq. (2.21) with respect to Lebesgue measure m suggests 


1 n —1 

- £</W> ra ~ </>„• • (2.22) 

i= 0 

Note that for z < 2 one has to choose a = 1, because the map still exhibits an SR.B measure, 
and Eq. (2.22) becomes an equality. However, for z > 2 we have an infinite invariant measure 
that cannot be normalized, hence here Eq. (2.22) remains a proportionality, unless we fix 
this constant by other constraints. 

These known facts from infinite ergodic theory motivate to suitably define generalized chaos 
quantities, which assess weakly chaotic dynamics by yielding non-zero values. Following 
the left hand side of Eq. (2.22), by choosing f(x) = In \P' a z {x)\ we define the generalized 
Lyapunov exponent as 

A := lim r(1 t a) E( hl l M/ (^)l)m • (2.23) 

n—>oo 77, a z ' 

i =0 

The inclusion of the gamma function in the numerator is not obvious at this point, however, 
it turns out to be convenient when calculating A for the Pomeau-Manneville map [How09]. 
Interestingly, it is precisely the same canonical choice as is made in other areas of anomalous 
dynamics [Kor07]. Analogously, we amend Eq. (2.12) to define the generalized KS entropy 
as 

H ks := lim ~ r ^ ^ a ^ T; p*(w)lnp*(w) . (2.24) 

n—XX) n z —' 

we{wp} 
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Figure 2.4: Illustration of the interplay between weak chaos, infinite measures, and anoma¬ 
lous dynamics in the Pomeau-Manneville map Eq. (2.13) shown by the thick red lines: 
Anomalous dynamics is indicated by the irregular behaviour of the two single trajectories. 
Weak chaos is exemplified by the divergence of the two trajectories starting at nearby ini¬ 
tial conditions. The light blue area in the lower left corner depicts the shape of the infinite 
invariant density, which diverges at the marginal fixed point of the map [Klall]. 

Both quantities can be calculated independently for the piecewise linearization of P a , z pro¬ 
posed in Ref. [Gas88] by applying the thermodynamic formalism [Rue78, Bec93] in com¬ 
bination with transfer operator methods [Pre92, Tas02, Tas04]. As a result, one obtains 
[How09] 

H ks = A , (2.25) 

which may be considered as a generalization of Pesin’s formula Eq. (2.11) to anomalous 
dynamics. Related generalizations of chaos quantities, and other versions of a generalized 
Pesin formula, have been discussed in Refs. [Kor09, KorlO, AkilO, Saal2]. We remark, 
however, that in the mathematical literature there is the well-known formula by Rokhlin 
[Kel98], which for the Pomeau-Manneville map reads [Tha83, ZweOO] 

h Kr = (\n\P' az (x)\)^ . (2.26) 

Here the left hand side holds for the so-called Krengel entropy. In case of finite invariant 
measures one can show that Hks — hKr , the right hand is the Lyapunov exponent defined 
via the ensemble average Eq. (2.6), and Rokhlin’s formula boils down to Pesin’s formula 
Eq. (2.11). For infinite invariant measures, one can show that hx r — Hks [How09]. Com¬ 
bining Rokhlin’s formula with the integrated form of the Aaronson-Darling-Kac theorem 
Eq. (2.22) by using f{x) = ln|P' 2 (x)|, exploiting the definition Eq. (2.23) for the gen¬ 
eralized Lyapunov exponent, and by fixing the constant of proportionality in Eq. (2.22) 
with respect to Lebesgue initial measure, one recovers Eq. (2.25). One may thus argue 
that, within this setting, Rokhlin’s formula is a generalization of Pesin’s formula for infinite 
measure-preserving transformations, and that Eq. (2.25) is an illustration of it, worked out 
for the example of the Pomeau-Manneville map [How09]. 

The main theoretical objects of discussion in this subsection are shown together in Fig. 2.4. 
This figure actually represents the logo of the conference about Weak chaos, infinite ergodic 
theory and anomalous dynamics that was referred to in the introduction [Klall], from which 
the title of this book chapter derives. 
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2.3 A generalized hierarchy of chaos 

We conclude this section by embedding the previous results into the more general context 
of irregular deterministic dynamics [Kla07]. There exist in fact further fundamental types 
of dynamics that are intermediate between strongly chaotic, in the sense of exponential 
sensitivity quantified by a positive Lyapunov exponent, and trivially being non-chaotic 
in terms of purely regular dynamics. These different types of irregular dynamics can be 
characterized by suggesting a classification of chaotic behavior based on the dispersion of 
initially infinitesimally close trajectories. 

We start from the general expression for the asymptotic growth of the displacement A (t) of 
two trajectories generated by dynamics in continuous time t in the form of [ las88] 

In A (t) ~ t v °(\nt) Ul , 0 < v Q , zq G R. . (2.27) 

If Uq — 1, q = 0 we recover the usual exponential dynamical instability of Eq. (2.3), 

A(t) ~ exp (At) , (2.28) 

representing Lyapunov chaos [Ott93], whose strength is well quantified by the maximal 
positive Lyapunov exponent A. As discussed before, if A(i) grows weaker than exponential, 
one speaks of weak chaos [ZasOl, Gal03, vB04, Art05]. The regime of Eq. (2.27) with 0 < 
Uq < 1 or uq = 1 and zq < 0 , which is typical for intermittent dynamics, was characterized 
as sporadic by Gaspard and Wang [Gas88] ; cf. Eq. (2.17) and our respective discussion, 
as well as further details of this dynamics as presented in the following section. Here 
the dynamical instability is either of stretched exponential type or it is exponential with 
logarithmic corrections, 

A(t) ~exp(i I '°(lnf) 1 ' 1 ) . (2.29) 

Eq. (2.27) with u 0 = 0 and zq = 1, on the other hand, yields purely algebraic dispersion, 

A (t) ~ G 2 , 0 < zq , (2.30) 

for which Zaslavsky and Edelman [Zas02, Zas03] suggested the term pseudochaos. 2 Note 
that algebraic dispersion with logarithmic corrections may also exist, 

A(t) ~t U2 {\nty 3 ,zq G M , (2.31) 

covering a slightly larger class of dynamical systems. A prominent class of dynamical sys¬ 
tems exhibiting algebraic dispersion are polygonal billiards; two examples are depicted in 
Fig. 2.5. They represent the special case of pseudochaotic dynamics with zq = 1 for which 
the dispersion is strictly linear in time, 


A(t) ~ t . (2.32) 

However, in contrast to Lyapunov chaos and our weakly chaotic generalizations, here the 
linear dispersion does not actually capture the essential mechanism leading to dynamical 
randomness. For example, according to this classification both free flights and polygonal 

2 We remark that in Refs. [Zas02, Zas03, Zas05] one finds several slightly different definitions of pseu¬ 

dochaos. Here we refer to the first one stated in Ref. [Zas02]. 





Figure 2.5: Two simple examples of two-dimensional polygonal billiards [Ivla07]: A particle 
with unit velocity moves inside the depicted geometric domains by scattering elastically 
with their boundaries, (a) shows a right triangular billiard [Art97], (b) the triangle channel, 
where a unit cell with triangular scatterers is spatially continued along the line [Li02] . 


billiards of genus one, which clearly exhibit regular dynamics, are also pseudochaotic. As 
is discussed for the example of rational billiards, e.g., in Ref. [ vla07], in polygonal billiards 
complicated topologies reflecting the existence of corners, which yield pseudohyperbolic fixed 
points and pseudointegrability, provide the source of non-trivial irregular dynamics. One 
is thus tempted to speak of topology-induced chaos 3 as a subclass of pseudochaos if there 
is linear dispersion on surfaces that are not integrable. Pseudointegrahle rational billiards 
then form a subclass of topology-induced chaotic dynamical systems. Surprisingly, systems 
with linear dispersion generating non-trivial dynamics due to complicated topological struc¬ 
tures may still exhibit ergodic and transport properties as they are usually associated with 
Lyapunov unstable chaotic dynamical systems. The trivial end point of this attempt of a 
generalized classification of chaotic dynamics on the basis of dispersion is simply the purely 
regular, or periodic, case of A(f) = const. 


3 Anomalous diffusion 


We now establish a cross-link between weakly chaotic dynamics as discussed in the previous 
section and the problem of deterministic diffusion. The main question we address is what 
type of diffusion arises if we suitably spatially extend a simple dynamical system exhibiting 
anomalous dynamics. We first set up our model, which can be considered as a purely 
deterministic, anomalous version of a random walk on the line, and introduce the concept of 
anomalous diffusion. We then outline continuous time random walk theory, a standard tool 
in the theory of stochastic processes to study anomalous diffusion. The results obtained from 
this theory, worked out for our model, are compared to results from computer simulations. 
We conclude this chapter by deriving on the basis of this theory a generalized, fractional 


3 This should not be confused with topological chaos as defined in Ref. [Gas98]. 
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diffusion equation that approximately reproduces the probability density function (PDF) of 
our model. 


3.1 A simple model generating anomalous diffusion 

A straightforward way to define a spatially extended dynamical system based on the Pomeau- 
Manneville map discussed in Section 2.2 is as follows: By using 


P a , z (x) = x + ax z , 0 < x < - 


(3.1) 


of Eq. (2.13) without the modulus, as well as reflection symmetry, 

Pa,z(~ X ) = - P a,z( X ) i (3.2) 

we continue this map onto the whole real line by a lift of degree one [Fuj82, Gei82, Sch82] 


Pa,z( X + !) = Pa,z{ X ) + 1 ' 


(3-3) 


The resulting model [Gei84, Zum93] is displayed in Fig. 3.1. Here points are not restricted 
anymore onto the unit interval. Due to the coupling between different unit cells by eliminat¬ 
ing the modulus, there are now ‘jumps’ possible from unit interval to unit interval. One may 
thus think of this dynamical system as a fully deterministic, anomalous version of a simple 
random walk on the line. A basic question is now which type of diffusion is generated by 
this model? As usual, the diffusive behavior is quantified by the mean square displacement 
(MSD) defined by 



where (x 2 ) is the second moment of the position PDF p n (x) at time step n. Starting from 
a given initial PDF po(x) at time step n = 0, points, or point particles, will spread out 
over the whole real line, as quantified by p n {x). Surprisingly, by calculating this MSD both 
analytically and from computer simulations one finds [Gci84, Zum93] that for z > 2 

(x 2 ) ~ n a , a < 1 (n —» oo) . (3.5) 


If one defines the diffusion coefficient of the system in the standard way by 


D := 



n—toc 277. ’ 


(3.6) 


Eq. (3.5) implies that D — 0, despite the fact that particles can go anywhere on the real line 
as illustrated in Fig. 3.1. While a process like Brownian motion leads to D > 0, here we thus 
encounter a very different type of diffusion: If the exponent a in the temporal spreading of 
the MSD Eq. (3.5) of an ensemble of particles is not equal to one, one speaks of anomalous 
diffusion [MetOO, Kla08]. If a < 1 one says that there is subdiffusion , for a > 1 there is 
superdiffusion. In case of linear spreading with a = 1 one refers to normal diffusion. The 
constant 


K ■ = 


, < x 2 > 

lim - 

n—Kx) Tl a 




(3.7) 
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Figure 3.1: The Pomeau-Mannevillc map Fig. 2.2, Eq. (2.13), lifted symmetrically onto the 
whole real line such that it generates subdiffusion. 

where in case of normal diffusion in one dimension K = 2 D, is called the generalized 
diffusion coefficient. 1 For our simple map model depicted in Fig. 3.1 we will first calculate 
the MSD analytically by means of stochstic continuous time random walk (CTRW) theory. 
By comparing the analytical formula with results from computer simulations, we will then 
focus on how K behaves as a function of z for fixed a revealing some interesting, non-trivial 
properties. 


3.2 Continuous time random walk theory 

Pioneered by Montroll, Weiss and Scher [Mon65, Mon73, Sch75], CTRW theory yields per¬ 
haps the most fundamental theoretical approach to explain anomalous diffusion [Bou90, 
Wei94, Ebe05]. In further groundbreaking works by Geisel et al. and Klafter et al., this 
method was then adapted to sub- and superdiffusive deterministic maps [Gei84, Gei85, 
Shl85, Zum93] 

The basic assumption of this approach is that diffusion can be decomposed into two stochas¬ 
tic processes characterized by waiting times and jumps, respectively. Thus one has two 
sequences of independent identically distributed random variables, namely a sequence of 


1 In detail, the definition of a generalized diffusion coefficient is a bit more subtle [Kor07]. 
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3 Anomalous diffusion 


positive random waiting times T\, T 2 , T 3 ,... with PDF w(t) and a sequence of random jumps 
Ci, C 2 , C 3 ? - - - with a PDF A(x). For example, if a particle starts at point x = 0 at time t 0 = 0 
and makes a jump of length ( n at time t n = Tf + T 2 + ... + T n , its position is x = 0 for 
0 < t < Tj = ti and x = Ci + C 2 + ••• + Cn for t n < t < t n+ 1 . The probability that at least 
one jump is performed within the time interval [0, t) is then f Q dt'w(t') while the probability 
for no jump during this time interval reads '&(£) = 1 — f* dt'w(t'). The master equation for 
the PDF P(x, t ) to find a particle at position x and time t is then 

/ OO pt 

dx'X(x — x') / dt' w(t — t' ) P(x > , t') + ^(t)5(x) , (3.8) 

-00 J 0 

which has the following probabilistic meaning: The PDF to find a particle at position x at 
time t is equal to the PDF to find it at point x' at some previous time t' multiplied with 
the transition probability to get from (x', t') to (x, t ) integrated over all possible values of 
x' and t'. The second term accounts for the probability of remaining at the initial position 
x = 0. The most convenient representation of this equation is in Fourier-Laplace space, 

^ poo poo 

P(k,s)= dxe ikx dt e~ st P(x,t) , (3.9) 

J —00 J 0 


where the hat stands for the Fourier transform and the tilde for the Laplace transform. 
Respectively, this function obeys the Fourier-Laplace transform of Eq. (3.8), which is called 
the Montroll-Weiss equation [Mon65, Mon73, Sch75], 


P(k,s) 


1 — w(s) 1 

s 1 — A (k)w(s) 


(3.10) 


The Laplace transform of the MSD can be obtained by differentiating the Fourier-Laplace 
transform of the PDF, 


<* 2 (*)> 



dx x 2 P(x, s ) 


d 2 P{k,s ) 
dk 2 


k =0 


(3.11) 


In order to calculate the MSD within this theory, it thus suffices to know A(x) and w(t) 
generating the stochastic process. For one-dimensional maps of the type of Eqs. (3.1),(3.2), 
by exploiting the symmetry of the map the waiting time distribution can be calculated from 
the approximation 


dx t z „ , 

Xn+1 - X n ~ — = ax z t , X < 1 


(3.12) 


where we have introduced the continuous time t > 0. This equation can be solved for x t 
with respect to an initial condition xq. Now one needs to define when a particle makes a 
“jump”, as will be discussed below. By inverting the solution for x t , one can then calculate 
the time t a particle has to wait before it makes a jump as a function of the initial condition 
Xq. This information determines the relation between the waiting time PDF w(t) and the 
as yet unknown PDF of injection points, 


dx 0 

dt 


m(t) — P in (x 0 ) 


(3.13) 
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Making the assumption that the PDF of injection points is uniform, Pj n ~ 1, the waiting 
time PDF is straightforwardly calculated from the knowledge of t(x o). The second ingredient 
that is needed for the CTRW approach is the jump PDF. Standard CTRW theory takes 
jumps between neighbouring cells only into account leading to the ansatz [Gei84, Zum93] 

A(x) = 5(|x| — 1) . (3-14) 

It turns out that in order to qualitiatively reproduce the dependence of the generalized 
diffusion coefficient K = K(z,a) Eq. (3.7) on the map’s two control parameters z and 
a, one needs to modify the standard theory at three points [Kor05, Kor07]: Firstly, the 
waiting time PDF must be calculated according to the unit interval [0,1], not according to 
[—0.5, 0.5], which is an alternative but not appropriate choice [Kla96b, Ivla97], yielding 

w(t) = a (1 + a(z — l)t) _S3T . (3.15) 


However, this PDF also accounts for attempted jumps to another cell, since after a step 
the particle may stay in the same cell with a probability of (1 — p). The latter quantity is 
roughly determined by the size of the escape region p — (1 — 2x c ) with x c as a solution of 
the equation x c + ax z = 1. We thus model this fact, secondly, by a jump length distribution 
in the form of 

M x ) = - 0 + (! ~P)8(x) . (3.16) 

Thirdly, in order to capture the dependence of K on z for fixed a, we define a typical jump 
length as 

i={|[M„(i)]|} , (3.17) 

where the square brackets stand for the floor function, which gives the coarse-grained dis¬ 
placement in units of elementary cells. The curly brackets denote both a time and ensemble 
average over points leaving a box. Note that for capturing the dependence of K on a for 
fixed z a different definition of the jump length is appropriate [Ivor05, Kor07]. Working 
out the modified CTRW approximation sketched above by taking these three details into 
account one obtains the result for the exponent a of the MSD, Eq. (3.5), 


a = 



\<z<2 


2 < z 


? 


(3.18) 


which matches to the standard theory [Gei84, Zum93]. This result is in excellent agreement 
to simulations for a broad range of control parameters z and a. Building on this result, the 
generalized diffusion coefficient can be calculated to 


K = pi 2 x 


a 7 sin (ny) / 7T7 1+7 
a(l - I/7) 


0 < 7 < 1 
1 < 7 < 00 


(3.19) 


with 7 := l/(z — l), which for z > 2 is identical to a of Eq. (3.18). In Fig. 3.2 this analytical 
approximation for K is compared with computer simulation results as a function of z for 
fixed a. There is good qualitative agreement between theory and simulations for z > 2, 
which converge asymptotically to each other for large z. For z < 2 there is reasonable 
qualitative agreement, though quantitative deviations, for z close to 2 while the theory does 
not work anymore for z —> 0, a problem that is discussed in Ref. [Kor07]. 
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Figure 3.2: The generalized diffusion coefficient K, Eq. (3.7), for the spatially extended 
Pomeau-Manneville map displayed in Fig. 3.1 as a function of z for a = 5. The bold 
black line depicts computer simulation results. The dashed line corresponds to the modified 
CTRW approximation Eqs. (3.17),(3. 19) [Kor07]. 


Remarkably, the K (z) obtained from simulations does not appear to be a smooth function of 
the control parameter, which is at variance with the prediction of CTRW theory Eq. (3.19). 
This non-smooth parameter dependence is not due to numerical errors (which here are 
very difficult to assess, as discussed below) but a well-known phenomenon for this type of 
systems. It has first been discovered for the normal diffusive case of this map with z = 1, 
where the diffusion coefficient D = K/2 has been studied both numerically and analytically 
as a function of the slope a as a control parameter [Kla95, Ivla99]. Note that for z — 1 the 
Pomeau-Manneville map boils down to a parameter-dependent, generalized version of the 
Bernoulli shift Eq. (2.2). 

We do not wish to further elaborate on the fractal parameter dependencies of transport 
coefficients in simple deterministic dynamical systems, an interesting phenomenon that has 
been discussed in detail in Refs. [Kla96b, Kla07]. Rather, we focus on the behaviour of 
the generalized diffusion coefficient at the point z — 2. According to the exponent a of 
the MSD given by Eq. (3.18), here the map exhibits a transition from normal to anomalous 
diffusion, which one may classify as a dynamical phase transition [Wan89, Bec93]. As can be 
seen in Fig. 3.2, right around this transition point there are significant deviations between 
CTRW theory and the simulation results. Most notably, at z = 2 the CTRW approximation 
forms a non-differentiable little wedge by predicting K{2) = 0, whereas the simulations yield 
K{2) > 0. By increasing the computation time one indeed finds very slow convergence of 
the simulation data towards the CTRW solution [Kor07]. 

The explanation of these deviations, and of the phenomenon of a complete suppression of 
the generalized diffusion coefficient right at the transition point, is obtained by carrying 
out a refined analysis by means of CTRW theory. For a long time it was known already 
that at z = 2, the MSD behaves like ( x 2 ) ~ n/lnn (n — Y oo) [Gei84, Zum93]. Note that 
according to our definition of the generalized diffusion coefficient Eq. (3.7) this logarithmic 
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depencence was incorporated into the strength of the diffusion coefficient, otherwise our 
analytical CTRW approximation would not have been continuous at z = 2. By taking 
into account higher-order terms when performing the CTRW theory calculations, which 
correspond to lower-order terms in time for the MSD, one arrives at 

, n < n cr and ~ n , n n cr , z < 2 

!% > ‘ = 2 ( 3 - 2 °) 

^ , n < n cr and ~ n“ , n > h cr , z > 2 . 

Here n cr and n cr are crossover times that can be calculated exactly in terms of the map’s 
control parameters. For z —» 2 both these crossover times diverge, and one arrives at the 
asymptotic nj Inn dependence mentioned before. The perhaps surprising fact is that around 
the transition point these multiplicative logarithmic corrections still survive for long but fi¬ 
nite time, in agreement with computer simulation results. In other words, these logarithmic 
corrections lead to an ultraslow convergence of the simulation results thus explaining the 
deviations between long-time CTRW theory and simulations shown in Fig. 3.2. But more 
importantly, these logarithmic terms dominate the strength of the generalized diffusion co¬ 
efficient around the transition point from normal to anomalous diffusion eventually yielding 
a full suppression of this quantity right at the transition point. It can be conjectured that 
the presence of such multiplicative logarithmic corrections around transition points between 
normal and anomalous diffusion is a typical scenario in this type of systems [Kor07]. 


3.3 A fractional diffusion equation 


We now turn to the PDFs generated by the lifted map Eq. (2.13). As we will show now, 
CTRW theory not only predicts the power a correctly but also the form of the coarse grained 
PDF P(x,t ) of displacements. Correspondingly the anomalous diffusion process generated 
by our model is not described by an ordinary diffusion equation but by a generalization of it. 
Starting from the Montroll-Weiss equation and making use of the expressions for the jump 
and waiting time distributions Eqs. (3.14), (3.15), we rewrite Eq. (3.10) in the long-time 
and -space asymptotic form 

sTP - s' 7 ” 1 = — ^—k 2 P (3.21) 

2 cff V ; 

with c = T(1 —7) and b = 7 /a. For the initial condition P(x, 0) = 5(x) of the PDF we have 

P(k, 0) = 1. Interestingly, the left hand side of this equation corresponds to the definition 

of the Caputo fractional derivative of a function G [Pod99, Mai97], 


<9 7 G 

dP 


in Laplace space, 


dt e 


r(i -7) 

CPG 


—st 


dP 


= s 7 G(s) - s 7 " 1 G(0) 


,dG 
dt' ’ 

(3.22) 

■G( 0) • 

(3.23) 


Thus, fractional derivatives come naturally into play as a suitable mathematical formalism 
whenever there are power law memory kernels in space and/or time generating anomalous 
dynamics; see, e.g., Refs. [Sok02, MetOO] for short introductions to fractional derivatives 








Figure 3.3: Comparison of the probability density obtained from simulations of the lifted 
map Eq. (2.13) (oscillatory structure) with the analytical solution Eq. (3.25) of the fractional 
diffusion equation Eq. (3.24) (continuous line in the middle) for z = 3 and a = 8 . The 
probability density was computed from 10' particles after n = 10 3 iterations. For the 
generalized diffusion coefficient in Eq. (3.25) the simulation result was used. The crosses 
(x) represent the numerical results coarse grained over unit intervals. The upper and the 
lower curves correspond to fits with a stretched exponential and a Gaussian distribution, 
respectively. The inset depicts the probability density function for the map on the unit 
interval with periodic boundaries. 


and Ref. [Pod99] for a detailed exposition. Turning back now to real space, we thus arrive 
at the time-fractional diffusion equation 


d 1 P(x, t) _ d 2 P 
dV dx 2 


(3.24) 


with D = AT(1 + y)/2, 0 < 7 < 1, which is an example of a fractional diffusion equation 
generating subdiffusion. Note the existence of the gamma function in the numerator defining 
D, which is analogous to the appearance of the gamma function in our generalized chaos 
quantities Eqs. (2.23),(2.24). For 7 = 1 we recover the ordinary diffusion equation. The 
solution of Eq. (3.24) can be expressed in terms of an M-function of Wright type [Mai97] 
and reads 

= 2y/Dtrft M (^2) ' (3 ‘ 25) 

Fig. 3.3 demonstrates an excellent agreement between the analytical solution Eq. (3.25) and 
the PDF obtained from simulations of the map Eq. (2.13) if it is coarse grained over unit 
intervals. However, it also shows that the coarse graining eliminates a periodic fine structure 
that is not captured by Eq. (3.25). This fine structure derives from the ‘microscopic’ invari¬ 
ant density of an elementary cell (with periodic boundaries) as shown in the inset of Fig. 3.3 
[Kla96b]. The singularities are due to the marginal fixed points of the map, where points 
are trapped for long times. Remarkably, that way the microscopic origin of the intermittent 
dynamics is reflected in the shape of the PDF on the whole real line: From Fig. 3.3 it is seen 
that the oscillations in the density are bounded by two functions, the upper curve being of 
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a stretched exponential type while the lower is Gaussian. These two envelopes correspond 
to the laminar and chaotic parts of the motion, respectively. 


4 Anomalous fluctuation relations 


After having accomplished a transition from deterministic dynamics to stochastic modeling 
in the previous section, for the reminder of this chapter we fully focus on stochastic sys¬ 
tems. First, we discuss a remarkable finding in nonequilibrium statistical mechanics that 
was widely investigated over the past two decades, which are fluctuation relations. After 
providing a brief outline of what they are, and why they are important, we first study an 
example of them for one of the most simple types of stochastic dynamics, which is Brown¬ 
ian motion modeled by an ordinary Langevin equation. Along these lines, we then consider 
generalized versions of Langevin dynamics exhibiting anomalous diffusion. For these types 
of dynamics we check again for fluctuation relations and in one case obtain a different, new 
form of such a formula. We argue that generalized, anomalous fluctuation relations should 
be important to understand nonequilibrium fluctuations in glassy dynamics. 


4.1 Fluctuation relations 

Fluctuation Relations (FRs) denote a set of symmetry relations describing large-deviation 
properties of the PDFs of statistical physical observables far from equilibrium. First forms 
defining one subset of them, often referred to as Fluctuation Theorems , emerged from gener¬ 
alizing fluctuation-dissipation relations to nonlinear stochastic processes [Boc81a, Boc81b]. 
They were then discovered as generalizations of the Second Law of Thermodynamics for ther- 
mostated dynamical systems, i.e., systems interacting with thermal reservoirs, in nonequilib¬ 
rium steady states [Eva93, Eva94, Gal95a, Gal95b]. Another subset, called work relations , 
generalize a relation between work and free energy, known from equilibrium thermodynam¬ 
ics, to nonequilibrium situations [Jar97b, Jar97a], These two fundamental classes were later 
on amended and generalized by a variety of other FRs from which they can partially be 
derived as special cases [Cro99, HatOl, Sei05, SaglO]. Research performed over the past 
ten years has shown that FRs hold for a great variety of systems thus featuring one of the 
rare statistical physical principles that is valid even very far from equilibrium; see, e.g., 
Refs. [Gal98, Eva02, Kla07, Har07, Sci08, Jar08, I\lal2b] and further references therein. 
Many of these relations have meanwhile been verified in experiments on small systems, i.e., 
systems on molecular scales featuring only a limited number of relevant degrees of freedom 
[Wan02, Rit03, Bus05, CillO, ToylO, Alell]. 

So far FRs have mostly been studied for dynamics exhibiting normal diffusion. This raises 
the question to which extent the ‘conventional’ FRs derived for these cases are valid for 
anomalous dynamics [Che09, Klal2a], Theoretical results for generalized Langevin equa¬ 
tions [Bcc04, Ohk07, Mai07, Cha08], Levy flights [Tou07, Tou09] and continuous-time ran¬ 
dom walk models [Esp08] as well as computer simulations for glassy dynamics [Sel09] showed 
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Figure 4.1: Sketch of a colloidal particle confined within a harmonic trap that is dragged 
through water with a constant velocity v*, cf. the experiment by Wang et al. [Wan02]. 


both validity and violations of the various types of conventional FRs referred to above, de¬ 
pending on the specific type of anomalous dynamics considered and the nonequilibrium 
conditions that have been applied [Che09]. 

In this section we outline how the two different fields of FRs and anomalous dynamics can 
be cross-linked in order to explore to which extent conventional forms of FRs are valid for 
anomalous dynamics. With the term anomalous fluctuation relations we thus refer to devia¬ 
tions from conventional forms of FRs, which are due to anomalous dynamics. Here we focus 
on generic types of stochastic anomalous dynamics by only checking transient fluctuation 
relations (TFRs), which describe the approach from a given initial distribution towards a 
(non)equilibrium steady state. As a warm-up, we first derive the conventional TFR for the 
trivial case of Brownian motion of a particle moving under a constant external force modeled 
by standard Langevin dynamics. We then consider a straightforward generalization of this 
type of dynamics in form of long-time correlated Gaussian stochastic processes. For two 
fundamental, different versions of this dynamics we check for the existence of conventional 
TFRs under the simple nonequilibrium condition of a constant external force. 

4.2 Fluctuation relations for ordinary Langevin dynamics 

Consider a particle system evolving from some initial state at time t — 0 into a nonequilib¬ 
rium steady state for t —> oo. A famous example that has been investigated experimentally 
[Wan02] is a colloidal particle immersed into water and confined by an optical harmonic 
trap, see Fig. 4.1. The trap is first at rest but then dragged through water with a constant 
velocity v*. 

The key for obtaining FRs in such systems is to compute the PDF p(£t) of suitably defined 
dimensionless entropy production over trajectory segments of time length t. The goal is 
to quantify the asymmetry between positive and negative entropy production in p(£t) f° r 
different times t since, as we will demonstrate in a moment, this relation is intimately related 
to the Second Law of Thermodynamics. For a very large class of systems, and under rather 
general conditions, it was shown that the following equation holds [Eva02, Har07, Kur07]: 

ln FW = ?i ' (41) 

Given that here we consider the transient evolution of a system from an initial into a steady 
state, this formula became known as the transient fluctuation relation (TFR). The left hand 
side we may call the fluctuation ratio. Relations exhibiting this functional form have first 
been proposed in the seminal work by Evans, Cohen and Morriss [Eva93], although in the 
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Figure 4.2: Illustration of the dynamics of the probability density function for entropy 
production p(£t) for different times t\ < t 2 < t 3 . 

different situation of considering nonequilibrium steady states. Such a steady state relation 
was proved a few years later on by Gallavotti and Cohen for deterministic dynamical systems, 
based on the chaotic hypothesis [Gal95a, Gal95b]. The idea to consider such relations for 
transient dynamics was first put forward by Evans and Searles [Eva94]. 

Fig. 4.2 displays the temporal evolution of the PDF for entropy production in such a situ¬ 
ation. The asymmetry of the evolving distribution, formalized by the fluctuation relation 
Eq. (4.1), is in line with the Second Law of Thermodynamics. This easily follows from 
Eq. (4.1) by noting that 

p(Zt) = p(~€t) exp(&) > p(-Ct) , (4.2) 

where is taken to be positive or zero. Integration from zero to infinity over both sides of 
this inequality after multiplication with and defining the ensemble average over the given 
PDF as (...) = J"*_ t/(, />/(,]... yields 


(6) > 0 . (4.3) 

As a warm-up, we may first check the TFR for the ordinary overdamped Langevin equation 
[Kub92] 

x = F + C(t) , (4.4) 

with a constant external force given by F and Gaussian white noise ((t). Note that for sake 
of simplicity, here we set all the other constants that are not relevant within this specific 
context equal to one. For Langevin dynamics with a constant force the entropy production 
defined by the heat, or equivalently the dissipative work, is simply equal to the mechanical 
work [vZ03b] 

W t = Fx(t) . (4.5) 

It follows that the PDF of entropy production, which here is identical to the one for the 
mechanical work, is trivially related to the PDF of the position x of the Langevin particle 
via 

p(W t )=F- 1 e(x,t) . (4.6) 

This is very convenient, since it implies that all that remains to be done in order to check 
the TFR Eq. (4.1) is to solve the Fokker-Planck equation for the position PDF g(x,t) f° r 
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a given initial condition. Here and in the following, we choose x(0) = 0, i.e., in terms of 
position PDFs we start with a delta-distribution at x — 0. Note that for ordinary Langevin 
dynamics in a given potential, typically the equilibrium density is taken as the initial density 
[vZ03a, vZ03b]. However, since in the following we will consider dynamics that may not 
exhibit a simple equilibrium state, without loss of generality here we make a different choice. 
For the ordinary Langevin dynamics Eq. (4.4) modeling a linear Gaussian stochastic process, 
the position PDF is Gaussian exhibiting normal diffusion [Kub92, Ris96], 


e(x,t) 



( ( *-{x)) 2 \ 

V 2<o ) 


(4.7) 


With the subscript zero we denote ensemble averages in case of zero external field. By using 
the PDF-scaling Eq. (4.6) and plugging this result into the TFR Eq. (4.1), we easily derive 
that the TFR for the work W t holds if 

(Wt) = % , (4-8) 

which is nothing else than an example of the fluctuation-dissipation relation of the first 
kind (FDR1) [Ivub66, Kub92]. We thus arrive at the seemingly trivial but nevertheless 
important result that for this simple Gaussian stochastic process, the validity of FDR1 
Eq. (4.8) implies the validity of the work TFR Eq. (4.1). For a full analysis of FRs of 
ordinary Langevin dynamics we refer to van Zon and Cohen Refs. [vZ03a, vZ03b]. 
Probably inspired by the experiment of Ref. [Wan02], typically Langevin dynamics in a 
harmonic potential moving with a constant velocity has been studied in the literature 
[Zam05a, Mai07, Ohk07, Cha08], cf. Fig. 4.1. Note that in this slightly more complicated 
case the (total) work is not equal to the heat [vZ03b]. While for the work one recovers 
the TFR in its conventional form Eq. (4.1) in analogy to the calculation above, surpris¬ 
ingly the TFR for heat looks different for large enough fluctuations. The origin of this 
phenomenon has been discussed in detail in Ref. [vZ03a], related effects have been reported 
in Refs. [Har06, Eva05, Har07]. However, in the following we check for another source of 
deviations from the conventional TFR Eq. (4.1) that is due to the existence of microscopic 
correlations in form of anomalous dynamics. In order to illustrate the main ideas it suffices 
to consider a nonequilibrium situation simply generated by a constant external force. 


4.3 Fluctuation relations for anomalous Langevin 
dynamics 

In our presentation of this section we follow Ref. [Che09], which may be consulted for further 
details. Our goal is to check the TFR Eq. (4.1) for Gaussian stochastic processes generating 
anomalous diffusion. These processes are defined by the overdamped generalized Langevin 
equation 

dt'x(t')^[t — t') — F + ((t) (4.9) 

with Gaussian noise ((t) and friction that is modeled with a memory kernel 7 (f). By us¬ 
ing this equation a stochastic process can be defined that exhibits normal statistics but 
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with anomalous memory properties in form of non-Markovian long-time correlated Gaus¬ 
sian noise. Equations of this type can be traced back at least to work by Mori and Kubo 
around 1965 (see [Ivub66] and further references therein). They form a class of stan¬ 
dard models generating anomalous diffusion that has been widely investigated, see, e.g., 
Refs. [Kub92, Por96, LutOl]. FRs for this type of dynamics have more recently been an¬ 
alyzed in Refs. [Bec04, Ohk07, Mai07, Cha08]. Examples of applications for this type 
of stochastic modeling are given by generalized elastic models [TallO], polymer dynamics 
[PanlO] and biological cell migration [Die08]. We split this class into two specific cases: 

1. Correlated internal Gaussian noise 

We speak of internal Gaussian noise in the sense that we require the system to exhibit the 
fluctuation-dissipation relation of the second kind (FDR2) [Kub66, Ivub92] 

mat)) ~ y(t -1), (4.io) 

again by neglecting all constants that are not relevant for the main point we wish to make 
here. We now consider the specific case that both the noise and the friction are correlated 
by a simple power law, 

y(t) ~ t -/3 , 0 < fl < 1. (4-11) 

Because of the linearity of the generalized Langevin equation (4.9) the position PDF must 
be the Gaussian Eq. (4.7), and by the scaling of Eq. (4.6) we have p{W t ) ~ g{x,t). It 
thus remains to solve Eq. (4.9) for mean and variance, which can be done in Laplace space 
[Che09] yielding subdiffusion , 

a x,F ~ ^ ■> (4-12) 

by preserving the FDR1 Eq. (4.8). Here and in the following we denote ensemble averages 
in case of a non-zero external held with the subscript F. For Gaussian stochastic processes 
we have seen in the previous section that the conventional work TFR follows from FDR1. 
Hence, for the above power-law correlated internal Gaussian noise we recover the conven¬ 
tional work TFR Eq. (4.1). 

2. Correlated external Gaussian noise 

As a second case, we consider the overdamped generalized Langevin equation 

x = F + C(t), (4.13) 

which represents a special case of Eq. (4.9) with a memory kernel modeled by a delta- 
function. Again we use correlated Gaussian noise defined by the power law 

(CWC(0> ~ I* -<1A o</?< 1, (4.14) 

which one may call external , because in this case we do not postulate the existence of 
FDR2. The position PDF is again Gaussian, and as before p(W t ) ~ g(x,t). However, by 
solving the Langevin equation along the same lines as in the previous case, here one obtains 
superdiffusion by breaking FDR1, 

t 2 ~P. 


fW t ) ~ t , <X\Vt,F 




(4.15) 
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Figure 4.3: The fluctuation ratio ln(II T (J T )/II T (J T )) for the entropy production W T = EJ r 
with particle current J T and field strength E for particle density p at different times r. 
The full line, with slope one, displays the result of the conventional FR Eq. (4.1) in a 
nonequilibrium steady state. The figure is from Ref. [Sel09]. 


Calculating the fluctuation ratio, i.e., the left hand side of Eq. (4.1), from these results 
yields the anomalous work TFR 

111 -T^k = C ^~ lw t 0 < p < 1, (4.16) 

where Cp is a constant that depends on physical parameters [Che09]. Comparing this equa¬ 
tion with the conventional form of the TFR Eq. (4.1) one observes that the fluctuation ratio 
is still linear in Wt thus exhibiting an exponential large deviation form [Tou09]. However, 
there are two important deviations: (1) the slope of the fluctuation ratio as a function of 
Wt is not equal to one anymore, and in particular (2) it decreases with time. We may thus 
classify Eq. (4.16) as a weak violation of the conventional TFR. 

We remark that for driven glassy systems FRs have already been obtained displaying slopes 
that are not equal to one. Within this context it has been suggested to capture these 
deviations from one by introducing the concept of an ‘effective temperature’ [Sel98, Zam05b, 
Zam05a]. As far as the time dependence of the coefficient is concerned, such behavior has 
recently been observed in computer simulations of a paradigmatic two dimensional lattice gas 
model generating glassy dynamics [Sel09] . Fig. 4.3 shows the fluctuation ratio as a function 
of the entropy production at different times r as extracted from computer simulations of 
this model, where the PDF has first been relaxed into a nonequilibrium steady state. It is 
clearly seen that the slope decreases with time, which is in line with the prediction of the 
anomalous TFR Eq. (4.16). To which extent the nonequilibrium dynamics of this lattice 
gas model can be mapped onto the generalized Langevin equation Eq. (4.13) is an open 
question. 

In summary, for Gaussian stochastic processes with correlated noise the existence of FDR2 
implies the existence of FDR1, and FDR1 in turn implies the existence of work TFR in 
conventional form. That is, the conventional work TFR holds for internal noise. However, 
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there is a weak violation of the conventional form in case of external noise yielding a pre¬ 
factor that is not equal to one and in particular depends on time. 


5 Anomalous dynamics of biological 
cell migration 


In order to illustrate the importance of anomalous dynamics for realistic situations, in this 
final section of our book chapter we discuss experiments and theory about the migration 
of single biological cells crawling on surfaces as an example. Here we focus on cells in an 
equilibrium situation, i.e., not moving under the influence of any external gradients or fields. 
This case is investigated by extracting results for the MSD and for the position PDFs from 
experimental data. We then show how the experimental results can be understood by a 
mathematical model in form of a fractional Klein-Kramers equation. As far as MSD and 
velocity autocorrelation function are concerned, this equation bears some similarity to a 
generalized Langevin equation that is of the same type as the one that has been discussed 
in Section 4.3. Our presentation in this section is based on Ref. [Die08]. 


5.1 Cell migration 

Nearly all cells in the human body are mobile at a given time during their life cycle. Em- 
bryogenesis, wound-healing, immune defense and the formation of tumor metastases are well 
known phenomena that rely on cell migration [Lau96, Lam09, FrilO], Fig. 5.1 depicts the 
path of a single biological cell crawling on a substrate measured in an in vitro experiment 
[Die08]. At first sight, the path looks like the trajectory of a Brownian particle generated, 
e.g., by the ordinary Langevin dynamics of Eq. (4.4). On the other hand, according to Ein¬ 
stein’s theory of Brownian motion a Brownian particle is passively driven by collisions from 
the surrounding fluid molecules, whereas biological cells move actively by themselves con¬ 
verting chemical into kinetic energy. This raises the question whether the random-looking 
paths of crawling biological cells can really be understood in terms of simple Brownian mo¬ 
tion [Dun87, Sto91] or whether more advanced concepts of dynamical modeling have to be 
applied [Har94, LIpaOl, Li08, Tak08, BodlO]. 


5.2 Experimental results and statistical analysis 

The cell migration experiments that we now discuss have been performed on two types of 
tumor-like migrating transformed renal epithelial Madin Darby canine kidney (MDCK-F) 
cell strains: wild-type ( NHE + ) and NHE-deficient ( NHE~) cells. Here NHE + stands 
for a molecular sodium hydrogen exchanger that either is present or deficient. It can thus 
be checked whether this microscopic exchanger has an influence on cell migration, which 



5 Anomalous dynamics of biological cell migration 



Figure 5.1: Overlay of a biological cell migrating in vitro on a substrate. The cell frequently 
changes its shape and direction during migration, as is shown by several cell contours ex¬ 
tracted during the migration process. The inset displays phase contrast images of the cell 
at the beginning and to the end of its migration process [Die08]. 


is a typical question asked by cell physiologists. The cell diameter is about 20-50/im and 
the mean velocity of the cells about lpm/min. Cells are driven by active protrusions of 
growing actin filaments ( lamellipodial dynamics ) and coordinated interactions with myosin 
motors and dynamically re-organizing cell-substrate contacts. The leading edge dynamics 
of a polarized cell proceeds at the order of seconds. Thirteen cells were observed for up to 
1000 minutes. Sequences of microscopic phase contrast images were taken and segmented 
to obtain the cell boundaries shown in Fig. 5.1; see Ref. [Die08] for full details of the 
experiments. 

According to the Langevin description of Brownian motion outlined in Section 4.2, Brownian 
motion is characterized by a MSD cr^ 0 (i) ~ t(t —* oo) designating normal diffusion. Fig. 5.2 
shows that both types of cells behave differently: First of all, MDCK-F NHE~ cells move 
less efficiently than NHE + cells resulting in a reduced MSD for all times. As is displayed 
in the upper part of this figure, the MSD of both cell types exhibits a crossover between 
three different dynamical regimes. These three phases can be best identified by extracting 
the time-dependent exponent /3 of the MSD cr^ 0 (t) ~ t 13 from the data, which can be done 
by using the logarithmic derivative 


m 


d In msd(t) 
dint 


(5.1) 


The results are shown in the lower part of Fig. 5.2. Phase I is characterized by an exponent 
/3(t) roughly below 1.8. In the subsequent intermediate phase II, the MSD reaches its 
strongest increase with a maximum exponent f3. When the cell has approximately moved 
beyond a square distance larger than its own mean square radius (indicated by arrows in 
the figure), j3(t) gradually decreases to about 1.4. Both cell types therefore do not exhibit 
normal diffusion, which would be characterized by j3(t) —» 1 in the long time limit, but move 
anomalously, where the exponent /3 > 1 indicates superdiffusion. 

We next study the PDF of cell positions. Since no correlations between x and y positions 
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Figure 5.2: Upper part: Double-logarithmic plot of the mean square displacement (MSD) 
as a function of time. Experimental data points for both cell types are shown by symbols. 
Different time scales are marked as phases I, II and III as discussed in the text. The solid 
lines represent fits to the MSD from the solution of our model, see Eq. (5.7). All parameter 
values of the model are given in [Die08]. The dashed lines indicate the uncertainties of the 
MSD values according to Bayes data analysis. Lower part: Logarithmic derivative /3{t) of 
the MSD for both cell types as defined by Eq. (5.1). 


could be found, it suffices to restrict ourselves to one dimension. Fig. 5.3 (a), (b) reveals 
the existence of non-Gaussian distributions at different times. The transition from a peaked 
distribution at short times to rather broad distributions at long times suggests again the 
existence of distinct dynamical processes acting on different time scales. The shape of these 
distributions can be quantified by calculating the kurtosis 


A) '■= 


(At)) 

(x 2 (t )) 2 


(5.2) 


which is displayed as a function of time in Fig. 5.3 (c). For both cell types K,(t) rapidly decays 
to a constant that is clearly below three in the long time limit. A value of three would be 
the result for the spreading Gaussian distributions characterizing Brownian motion. These 
findings are another strong manifestation of the anomalous nature of cell migration. 


5.3 Stochastic modeling 

We now present the stochastic model that we have used to reproduce the experimental data 
yielding the fit functions shown in the previous two figures. The model is defined by the 
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Figure 5.3: Spatio-temporal probability distributions P(x,t). (a),(b): Experimental data 
for both cell types at different times in semilogarithmic representation. The dark lines, 
labeled FKK, show the long-time asymptotic solutions of our model Eq. (5.3) with the same 
parameter set used for the MSD fit. The light lines, labeled OU, depict fits by the Gaussian 
distributions Eq. (4.7) representing Brownian motion. For t = 1 min both P(x,t ) show a 
peaked structure clearly deviating from a Gaussian form, (c) The kurtosis n(t) of P(x,t), 
cf. Eq. (5.2), plotted as a function of time saturates at a value different from the one of 
Brownian motion (line at k — 3). The other two lines represent n{t) obtained from the 
model Eq. (5.3) [Die08]. 


fractional Klein-Kramers equation [BarOO] 


dg d r , d 1 “ 
dt = ^ + dt^ la 
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(5.3) 


Here g = g(x,v,t ) is the PDF depending on time t, position x and velocity v in one 
dimension, y Q is a friction term and v^ h = ksT /M stands for the thermal velocity squared of 
a particle of mass M = 1 at temperature T, where ks is Boltzmann’s constant. The last term 
in this equation models diffusion in velocity space. In contrast to Fokker-Planck equations, 
this equation features time evolution both in position and velocity space. What distinguishes 
this equation from an ordinary Klein-Kramers equation, the most general model of Brownian 
motion [Ris96], is the presence of the Riemann-Liouville fractional derivative of order 1 — a 
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in front of the terms in square brackets. Note that for a = 1 the ordinary Klein-Kramers 
equation is recovered. The analytical solution of this equation for the MSD has been calcu¬ 
lated in Ref. [BarOO] to 

<7x,o(*) = 2^f 2 £ a>3 (-7 a f a ) -> 2 Y{3- a ) ^ °°) ( 5 - 5 ) 

with D a = v 2 h /^ a and the two-parametric or generalized Mittag-Leffler function (see, e.g., 
Chapter 4 of Ref. [Kla08] and Refs. [Gor97, Pod99]) 

oo k 

EaJS< d = S ViZkTW) • a ’ P > °’ z €<c ■ < 5 '6) 

Note that E 1:1 (z) = exp ( 2 ), hence E aj/ g(z) is a generalized exponential function. We see 
that for long times Eq. (5.5) yields a power law, which reduces to the long-time Brownian 
motion result in case of a = 1 . 

In view of the experimental data shown in Fig. 5.2, Eq. (5.5) was amended by including the 
impact of random perturbations acting on very short time scales for which we take Gaussian 
white noise of variance r/ 2 . This leads to [Mar02] 

^l,o,noised) = . (5.7) 

The second term mimicks both measurement errors and fluctuations of the cell cytoskeleton. 
In case of the experiments with MDCK-F cells [Die08], the value of r) can be extracted from 
the experimental data and is larger than the estimated measurement error. Hence, this 
noise must largely be of a biological nature and may be understood as being generated by 
microscopic fluctuations of the lamellipodia in the experiment. 

The analytical solution of Eq. (5.3) for g(x,v,t) is not known, however, for large friction 7 ^ 
this equation boils down to a fractional diffusion equation for which g(x, t) can be calculated 
in terms of a Fox function [Sch89a]. The experimental data in Figs. 5.2 and 5.3 was then 
fitted consistently by using the above solutions with the four parameters Vth , a, 7 and rj 2 in 
Bayesian data analysis [Die08]. 

In summary, by statistical analysis of experimental data we have shown that the equilibrium 
migration of the biological cells under consideration is anomalous. Related anomalies have 
also been observed for other types of migrating cells [Har94, UpaOl, Li08, Tak08, BodlO]. 
Our experimental results are coherently reproduced by a mathematical model in form of 
a stochastic fractional equation. We now elaborate on possible physical and biological 
interpretations of our findings. 

First of all, we remark that the solutions of Eq. (5.3) for both the MSD and the velocity au¬ 
tocorrelation function match precisely to the solutions of the generalized Langevin equation 
[LutOl] 

v = — ( dt' 7 (t — t')v{t') + f(t) . (5.8) 

Jo 

Here f(t) holds for Gaussian white noise and 7 (t) ~ t~ a for a time-dependent friction co¬ 
efficient with a power law memory kernel, which alternatively could be written by using a 
fractional derivative [LutOl]. For 7 (f) ~ S(t) the ordinary Langevin equation is recovered. 
Note that the position PDF generated by this equation is Gaussian in the long time limit 





and thus does not match to the one of the fractional Klein-Kramers equation Eq. (5.3). How¬ 
ever, alternatively one could sample from a non-Gaussian £(i) to generate a non-Gaussian 
position PDF. Strictly speaking, despite equivalent MSD and velocity correlations Eqs. (5.3) 
and (5.8) define different classes of anomalous stochastic processes. The precise cross-links 
between the Langevin description and the fractional Klein-Kramers equation are subtle 
[EulOT] and to some extent still unknown. The advantage of Eq. (5.8) is that it allows 
more straightforwardly a possible biophysical interpretation of the origin of the observed 
anomalous MSD and velocity correlations, at least partially, in terms of the existence of a 
memory-dependent friction coefficient. The latter, in turn, might be explained by anoma¬ 
lous rheological properties of the cell cytoskeleton, which consists of a complex biopolymer 
gel [SemOT]. 

Secondly, what could be the possible biological significance of the observed anomalous cell 
migration? There is an ongoing debate about whether biological organisms such as, e.g., 
albatrosses, marine predators and fruit flies have managed to mimimize the search time 
for food in a way that matches to optimizing search strategies in terms of stochastic pro¬ 
cesses; see Refs. [Benll, Visll] and further references therein. In particular, it has been 
argued that Levy flights are superior to Brownian motion in order to fold sparsely, randomly 
distributed, replenishing food sources [Benll]. However, it was also shown that in other 
situations intermittent dynamics is more efficient than pure Levy motion [Benll], For our 
cell experiment, both the experimental data and the theoretical modeling suggest that there 
exists a slow diffusion on short time scales, whereas the long-time motion is much faster, 
which resembles intermittency as discussed in Ref. [Benll], Hence, the results on anomalous 
cell migration presented above might be biologically relevant in view of suitably optimized 
foraging strategies. 


6 Summary 


This chapter highlighted some fundamental aspects of anomalous dynamics: The scene was 
set by section 2, which reviewed basic ideas of weak chaos by establishing crosslinks to 
infinite ergodic theory. This branch of ergodic theory provides a rigorous mathematical 
approach to study weakly chaotic dynamical systems. In particular, we proposed suitable 
definitions of generalized chaos quantities assessing weakly chaotic dynamics by yielding a 
generalized version of Pesin’s theorem. We also outlined a generalized hierarchy of chaos 
on the basis of different functional forms of the dispersion exhibited by nearby trajectories 
of a deterministic dynamical system. In section 3 we related these concepts to the problem 
of anomalous diffusion by spatially extending our previously discussed simple map model. 
Applying stochastic continuous time random walk theory to this model in comparison to 
computer simulations, we learned about an intricate dynamical phase transition between 
normal and anomalous diffusion, governed by multiplicative logarithmic corrections in the 
mean square displacement. We also derived a fractional diffusion equation that reproduced 
the subdiffusive diffusive dynamics of this model on coarse scales. The subsequent sec¬ 
tion 4 elaborated on fluctuation relations, here understood as a large-deviation symmetry 
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property of the work probability distributions generated by a given stochastic dynamics far 
from equilibrium. We familiarized ourselves with the conventional form of transient work 
fluctuation relations derived from standard Langevin dynamics before exploring anomalous 
generalizations of Langevin equations. One of them reproduced the conventional form of 
fluctuation relations, whereas the other one yielded a generalized, anomalous fluctuation 
relation. The precise form of the resulting fluctuation relation appeared to be intimately 
releated to whether or not fluctuation-dissipation relations are broken. In our final main 
section we related our previous theoretical ideas to the experimental problem of studying 
biological cell migration. By extracting the mean square displacement and the position 
probability distributions from experimental data, we found that the dynamics exhibited by 
these cells was anomalous, showing different behavior on different time scales, by eventually 
yielding superdiffusion for long times. On the basis of these experimental results we sug¬ 
gested a stochastic theoretical model of cell migration in form of a fractional Klein-Kramers 
equation, which coherently reproduced our experimental findings. 

In summary, we traversed quite an anomalous scientific landscape of different but related 
topics: Starting from simple deterministic maps and their ergodic theory description we 
switched to basics of anomalous stochastic processes, studied both normal and anomalous 
stochastic fluctuations very far from equilibrium in terms of Langevin dynamics by ending 
up with anomalously crawling biological cells. We thus meant to illustrate the third column 
displayed in the very first Fig. 1.1 of the introduction, by also explaining the title of this 
contribution. Within a larger scientific context, one may consider our disucssion as an in¬ 
dication that a novel theory of anomalous nonequilibrium processes is presently emerging. 
In contrast to standard nonequilibrium statistical mechanics, this dynamics is inherently 
non-stationary, due to the weak chaos by which it is generated. This mechanism leads to 
important physical consequences like anomalous transport, which can be tested in experi¬ 
ments. On the side of theoretical physics this approach asks for further generalizations of 
recently developed fundamental concepts, perhaps leading to a weakly chaotic hypothesis, 
the identification of the physically relevant measures characterizing such systems, and to 
deriving experimentally measurable consequences such as generalizations of ordinary large 
deviation properties and fluctuation relations. However, these questions also motivate fur¬ 
ther mathematical work in upcoming directions of infinite ergodic theory to provide a formal 
framework and rigorous results for parts of the physical theory. 
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